Effective Rheology of Bubbles Moving in a Capillary Tube 
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We calculate the average volumetric flux versus pressure drop of bubbles moving in a single 
capillary tube with varying diameter, finding a square-root relation from mapping the flow equations 
onto that of a driven overdamped pendulum. The calculation is based on a derivation of the 
equation of motion of a bubble train from considering the capillary forces and the entropy production 
associated with the viscous flow. We also calculate the configurational probability of the positions 
of the bubbles. 

PACS numbers: 47.56.+r, 47.55.Ca, 47.55.dd, 89.75.Fb 



Multiphase flow in porous media plays a pivotal role 
in a vast range of applications in different fields such as 
oil recovery, soil mechanics and hydrology [D-Q • I n spite 
of its relevance in these important fields, fundamental 
questions still linger on. In particular, this is true in con- 
nection with steady-state multiphase flow which 
sets in after the initial instabilities such as viscous fin- 
gering are over in e.g. flooding experiments. A way to 
study this flow in the laboratory, is to simultaneously 
inject the two immiscible fluids into the porous medium 
and let them mix until a steady state where clusters and 
bubbles break up and mer ge b ut in such a way that their 
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averages remain constant 

Recently, the relation between average volumetric flow 
and excess pressure drop across the system has been in- 
vestigated both experimentally and theoretically fiol- 13|. 
The conclusion from these studies is that the volumetric 
flux depends quadratically on the excess pressure. This 
is in contrast to the assumptions of linearity commonly 
made when considering such systems, e.g. in connection 
with invoking the concepts of relative permeability in 
reservoir simulations at the flow rates where capillary 
and viscous forces compete 
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One aim of this work is to derive the volumetric flux 
versus excess pressure drop for a single capillary tube. 
We find that there is a square root singularity in this re- 
lation. This is in contrast to a the situation for network of 
pores, i.e., a porous medium. We base this calculation on 
the equation of motion of a bubble train in a long capil- 
lary tube with varying diameter, the Washburn equation 
[151 ] . We derive this equation following a different route 
than the now 91 year old original derivation. Lastly, we 
derive the probability distribution of the configuration of 
bubbles in the tube. This makes it possible to calculate 
the average of any quantity associated with the flow. 

We assume a long tube of length L oriented along an x 
axis. The radius of the tube, r, varies with the position 
along the tube as 
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FIG. 1: Shape of tube and corresponding pressure drops. 



where ro is the average radius of the tube, I is the wave- 
length along the tube and a is the dimcnsionless ampli- 
tude of the oscillation, sec Fig. [1] We assume L 3> /. 

We now imagine a bubble in a tube segment with 
length I limited by interfaces at Xi and Xf. The fluid 
in the bubble ("oil") is less wetting with respect to the 
tube walls than the fluid outside the bubble ("water"). 

The capillary pressure drop across the interface be- 
tween the two fluids at xa is [2fl 



2cr„ 



2<T„ 



1 



/ 2irx t 



r{xi) r 
and across the interface at Xf, 
2o~ wo < 2o~ wo 



r(xf) 



'o 



1 — a cos 



2ixxi 



(2) 



(3) 



Here, o~ wo is the surface tension between the two fluids. 
We sum the two capillary pressure drops to get 
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where we have defined Xb = (xi + Xf)/2 and Axb = 
xj — Xi . We assume the length of the bubble Axb to be 
smaller than the wavelength I, which is also equal to the 
length of the tube segment. Furthermore we have chosen 
the origin of the ^-coordinate such the whole bubble is 
located between and I. 

We assume the viscosities of the two fluids to be /j, w 
("water") and fi ("oil") respectively. In order to derive 
a constitutive equation between volumetric flux q and 
pressure drop along the tube, AP = (L/l)Ap, we con- 
sider the entropy production of a tube segment with an 
average cross-sectional area nr 2 , and length I [l6| . 
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where s(x) the entropy production in the tube per unit 
of length at position x. We assume the temperature to 
be constant along the tube. The entropy production 
in the fluid sections times the temperature is equal to 
—qdp(x)/dx, where q is the volumetric flux and p(x) the 
pressure at x. Because of the incompressible nature of 
the flow, the volumetric flux is independent of position. 
There is no excess entropy production at the interfaces 
between the fluid and the bubble, assuming equilibrium 
between the two immiscible fluids, cf. Eqs. @ and 
Substituting s(x) into Eq. ([5J and integrating gives 



dS 

= -q(Piw-pi)-q(pfo-Pio)-q{p2-Pfw) 



-Q (Ap - p c {x b )) 



(6) 
(7) 



where is the temperature, Ap = p-2 — p\ is the pressure 
drop across a length / and p c (xb) = Pio — Piw + Pfw — 
Pfo is the capillary pressure difference given in Eq. (j4|. 
The resulting linear force-flux relations for the pressure 
differences in the three sections arc written in the form 



Piw - Pi = 
Pfo - Pio = 

Pi - Pfw = 
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where R w and R Q are the Onsager resistivities per unit 
of length for the flow of water and bubble phases, respec- 
tively. Using Poisseulle flow the Onsager resistivities are 
equal to 8^/(7rrg), where \i is the viscosity of the fluid 
considered. The resulting pressure differences are 
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The sum of the pressure differences gives the constitutive 
equation for motion of a single bubble 
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where 
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is the average viscosity of the two fluids, which we note 
is independent of the position of the bubble. One may 
obtain a similar expression directly from Eq. (I7|. Eqs. 
(|T0| and ([TT]) were first derived by Washburn [15|, but 
by a different route. The present procedure clarifies why 
Eq. (flO|) should contain the average viscosity and how 
the contribution due to the capillary pressure arises. 

The center of mass coordinate of the bubble moves as 
irrgib = q. Hence, we have the equation of motion 
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We introduce the angle variable 6 = 2irxb/l and the time 
scale t = r ytTrrQ/(Al 2 p av ). We assume the pressure drop 
to be negative, Ap = — \Ap\. The equation of motion 
(|T3l then becomes 
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which is nothing but the equation of motion for the over- 
damped driven pendulum [171 ]. 

The period T T (measured in the same units as r), 
needed for the bubble to move from one end of the tube 
with length I to the other end for a given Ap, is 
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In seconds this time is 
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We then define the time-averaged angular speed as 
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Now, transforming back from dO/dr to ib to volumetric 
flux, q, we find the time-averaged flux equation 
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where sgn is the sign function and O is the Hcaviside 
function. Hence, for |Ap| < 7 there is no flow through 
the tube. Furthermore, if there is flow, T(xb) = I, as one 
would expect. 

We see that we can write Darcy's law for the time 
averaged volume flux only if 7 = 0, which is the case if 
the tube has a constant radius (a = 0). The deviation 
from this law for the time average is due to a capillary 
pressure which varies as a function of the position of the 
bubble along the tube. The effective flux equation (|18|) 
has a threshold pressure that must be overcome to induce 
flow, 7 defined in Eq. (|T3|) . Close to this threshold, when 
|Ap| - 7<7, the average flow equation becomes 



(?) = --^ 1 ^agn(Ap)e(\Ap\- 1 )^\A P \ - 1 , (19) 

i.e., there is a square root singularity. As shown in Rcf. 
17 1, this square root is a consequence of the quadratic 
cxtremum of r in Eq. (JTJ) leading to a saddlc-nodc bifur- 
cation, and therefore it does not depend on the specific 
sinusoidal shape of the profile chosen. The smaller the 
value of ro, the larger is the threshold value. A radius 
of 10 micrometer can give a threshold of the order of 80 
bar, which is non-negligible for most practical purposes. 

We may generalize these considerations to a tube seg- 
ment with length L — (27V + 1) Z in which there are 
27V + 1 bubbles numbered from —TV to TV. Each bub- 
ble j may be characterized by a center of mass position 
Xj and a width Axj. In view of the incompressible na- 
ture of the flow the volume flux q is independent of the 
position. This implies that the velocity Xj of all the bub- 
bles is the same, Xq = Xj. The equation of motion (fT2")) 
may then be generalized as follows, 
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where AP is the total pressure drop over the whole tube 
segment and 5xj = Xj — xq. Furthermore we define 
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By using trivial trigonometric identities, we may rewrite 
this expression as 
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It should be noted that T s and r c are proportional to the 
number of segments 27V + 1 with length I and therefore 
to the total length L. 

On non-dimensional form, Eq. (|22[) becomes 
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Choosing this form, we have assumed T s > T c . Hence, 
the saddle-node bifurcation will occur in the sine term 
and not in the cosine term, which may be ignored. Work- 
ing through the arguments leading to (fT8|) . we end up 
with the same effective flux equation as for the one- 
bubble case, but with 7 substituted for T s . If, on the 
other hand, T c > T s , we may shift 6 by tt/2, and we 
are back to Eq. (|25|) . but with T s and r c interchanged. 
Hence, again we find an effective flux equation as (|18|) . 
but with 7 substituted for r c . 

Hansen and Ramstad (lil ] proposed to approach 
steady-state immiscible two-phase flow in porous media 
using the methods of statistical mechanics. A central 
quantity in this context is the configurational probability 
njconfiguration}. Possessing this quantity makes it pos- 
sible to calculate the average of any quantity associated 
with the flow. We now derive the configurational prob- 
ability for bubbles in a capillary tube. The derivation is 
based on mapping the time average on to a configura- 
tional average. 

Time averaging assumes that the states in each time 
interval are equally probable. The state of the tube at 
time t is characterised by the position Xbit) of the bubble. 
The time average of a function f(xb(t)) is therefore given 
by 
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Using the relation between Xb and t we may write this 
average as an average over x b in the following way: 



7 V A P 2 "7 2 Jo jAr l +T^nlWi) dx b 



This may be written as 
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is the probability that the bubble has the position x b . 
This probability distribution can be interpreted as the 
probability distribution of an ensemble of tubes. Hence, 
it is the configurational probability. 
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The configurational probability depends on the man- 
ner the flow is controlled. If q is kept constant H(x b ) = 
(Tib)^ 1 = If Ap is kept constant one finds the value 
given on the right hand side of Eq. (j2U)) . Whether q or 
AP is kept constant is comparable to the choice of an 
ensemble. 

It is interesting to calculate a few averages in the con- 
stant AP ensemble. For the average velocity we find 
using Eqs.flU]) and j29j that 
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We calculate the average potential energy stored associ- 
ated with the capillary forces, using Eqs. (jTU)) . (f2l))) and 
Pc(x b ) = jsin(2irx b /l), finding 



(Pcq) = / U(x b )p c (x b )q(x b )dx b 
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The average potential energy associated with the capil- 
lary forces is therefore zero as it must be. 

We proceed to consider an ensemble of single tube seg- 
ments. For the ensemble we may interprete H(x b ) as the 
probability that the bubble has the position x b . For the 
ensemble of tubes this contributes fee \nlU(x b ) to the en- 
tropy density per unit of length along the tube. Together 
with the other entropy contributions in the single tube 
we then have 
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From thermodynamics it follows that the entropy den- 
sities per unit of volume of the single component bulk 
phases are given by s Q = —dp Q /dQ and s w = —dp w /dO. 
The reference value S° is due to other contributions to 
the entropy, like the excess entropies of the surfaces, 
which arc constant and which we do not need to consider 
explicitly. The last term, which is the configurational 
entropy, is constant. 

The configurational probability Eq. (|29p may be gen- 
eralized to a network of tubes, i.e., a porous medium [19j. 
This makes it possible to fulfill the program sketched by 
Hansen and Ramstad fl8j . 

We have in this note demonstrated that the average 
flux-pressure relation in a tube with varying diameter 



is non-linear and controlled by a square-root singularity. 
We have shown this by deriving the equation of motion 
from considering the capillary forces and the entropy pro- 
duction associated with the viscous flow. We have also 
derived the configurational probability, i.e., the probabil- 
ity of bubble configurations from the equation of motion. 
From this probability the average of any quantity associ- 
ated with the flow may be calculated. 

S. S. thanks the Norwegian Research Council for finan- 
cial support. 



[1 
[2 

[3; 
[4 
[5 
[6 

[r 
[8 
[9 
[10 



(32) [11 



[12 

[1.3 
[14 

[15 
[16 



[17 
[18 
[10 



* Electronic address: Santanu.Sinha@ntnu.no 
1 Electronic address: Alex.Hansen@ntnu.no 

* Electronic address: Dick.Bedeaux@ntnu.no 
Electronic address: Signe.Kjelstrup@ntnu.no 

] F. A. L Dullien, Porous Media — Fluid Transport and 

Pore Structure (Academic, New York, 1979). 
] J. Bear, Dynamics of Fluids in Porous Media (Dover 

Publ. Comp., Mineola, 1988). 
] P. M. Adler, Porous Media: Geometry and Transport 

(Butterworth-Heinemann, Stoneham, 1992). 
] M. Sahimi, Flow and Transport in Porous Media (VCH, 
Boston, 1995). 

] A. K. Gunstensen and D. H. Rothman, J. Geophys. Res. 

98, 6431 (1993). 
1] D. G. Avraam and A. C. Payatakes, J. Fluid Mech. 293, 
207 (1995); D. G. Avraam and A. C. Payatakes, Transp. 
Por. Media 20, 135 (1995). 
] D. G. Avraam and A. C. Payatakes, Ind. Eng. Chem. 

Res. 38, 778 (1999). 
1] H. A. Knudsen, E. Aker and A. Hansen, Transp. Por. 

Media, 47, 99 (2002). 
'] H. A. Knudsen and A. Hansen, Phys. Rev. E 65, 056310 
(2002). 

'] K. T. Tallakstad, H. A. Knudsen, T. Ramstad, G. L0voll, 
K. J. Mal0y, R. Toussaint and E. G. Flekk0y, Phys. Rev. 
Lett. 102, 074502 (2009). 
] K. T. Tallakstad, G. L0voll, H. A. Knudsen, T. Ramstad, 
E. G. Flekk0y and K. J. Mal0y, Phys. Rev. E, 80, 036308 
(2009). 

] E. M. Rassi, S. L. Codd and J. D. Seymour, New J. Phys. 

13, 015007 (2011). 
■J S. Sinha and A. Hansen, Europhys. Lett, in press (2012). 
] M. Carlson, Practical Reservoir Simulation (PennWell 

Corporation, Tulsa, 2006). 
] E. W. Washburn, Phys. Rev. 17, 273 (1921). 
1] S. Kjelstrup, D. Bedeaux, E. Johannessen and J. Gross, 
Non-Equilibrium Thermodynamics for Engineers (World 
Scientific, Singapore, 2010). 
] S. H. Strogatz, Non-Linear Dynamics and Chaos (Perseus 

Press, Cambridge, 1994). 
1] A. Hansen and T. Ramstad, Comp. Geosci. 13, 227 
(2009). 

'] S. Sinha, A. Hansen, D. Bedeaux and S. Kjelstrup, in 
preparation (2012). 



